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Dynamics of spin-glasses subjected to slow continuous changes of working enviroment 
such as slow changes of temperature or interaction bonds are studied based on scaling argu- 
ments and numerical simulations of continuous bond changes. Such perturbations lead to a 
glassy stationary state where the age or domain size of the system is pinned to macroscopic 
but finite values due to competition between relaxation and chaos effects (rejuvenation). 
Flutuation-dissipation relation is also pinned to that of a finite age. The scenario explains 
the anomalously weak cooling rate dependece of spin glasses. 

Introduction- Dynamics of glassy systems have been actively studied over the 
past several decades. One general interest has been the aging effect, 1 ) i.e. a glassy 
system ages extremely slowly in a fixed working environment specified for example 
by the temperature of the heatbath. Recently responses of such glassy systems to 
various kinds of perturbations have attracted much attention including restart of 
aging or rejuvenation of spin-glasses after sudden changes of the temperature of the 
heatbath, 2 )' 3 ^ responses of soft glassy systems to shear. 4 ) 

In the present work we study relaxation of spin-glasses subjected to slow changes 
of the working environment, namely slow changes of temperature or interaction 
bonds. This study is motivated by the recent experimental observations that spin- 
glasses exhibit surprisingly weak dependence on the cooling rates. 2 )' 3 ) We propose 
a scaling ansatz for the effective age of the system under continuous changes of the 
working environment assuming that the competition between aging and rejuvenation 
due to the chaos effects 5 ) leads to a stationary state. We performed Monte Carlo 
(MC) simulations of a 4-dimensional Edwards-Anderson (EA) model subjected to 
continuous changes of bonds and found good agreement with the scenario. 

Model- Specifically we consider the Edwards-Anderson (EA) spin-glass model 
with Ising spins S{ (i = 1. . . N) on a d-dimensional (hyper)-cubic lattice with N = L d 
lattice points. The Hamiltonian is given by H = — Yl<% - y JijSiSj — h ]TV Si. We use 
the ±J model where the interaction bonds J,j between nearest-neighbours take 
random dbJ values. The temperature is measured as k^T/J with the Boltzmann's 
constant set as &b = 1. The magnetic field h is measured as h/J. The relaxational 
dynamics is modeled by the usual heatbath single spin flip MC algorithm. 

A continuous bond changes is realized by flipping the sign of a fraction p C 1 
of the bonds chosen randomly in 1 Monte Cairo step (MCS). This amounts to flip a 
fraction p e s(t) = [1 — (1 — 2p)*]/2 of the original set of bonds Jij(0) to create the set 
of bonds Jij(t) at t (MCS). For pt < 1 p eS (t) = pt and p eS (t) 1/2 for pt » 1. 

A continuous temperature changes is specified by a cooling or heating rate vt = 
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6T/St where ST is the change of the temperature within an interval of time 5T. Thus 
the resultant change of the temperature after a time t is AT(t) = vrt. Unfortunately 
the temperature-chaos effect is much weaker than the bond chaos effect 6 -* in the 
present model and we limit ourselves to the bond perturbations in the simulations. 

Scaling arguments- Let us consider a spin-glass whose spins are evolving in 
a fixed working environment W = (J,T) specified by a set of bonds J = {Jij} 
and the heatbath temperature T. If the initial spin configuration is far from the 
equilibrium state associated with W, the system is equilibrated only up to a certain 
length scale Lx(t) after a time t. The domain size Lx(t) grows slowly with time t by 
thermally activated droplets excitations as L^(t) = Lq[(T/J) log(t/r c (T))] 1 /^ where 
Lq is a unit of length scale, r c (T) is the attempt time of the activated dynamics and 
if> > is the exponent of the free-energy barriers (See 7 ) for the details). 

Now we perturb the system by slowly changing W either by continuous tempera- 
ture changes characterized by AT(t) = vrt or continuous bond changes characterized 
by AJ(t) = J \Zpcs(t)-^ Within an interval of time r, the temporal working environ- 
ment W(t) changes from W(0) to W(t). Naturally the overlap length 5 ) between the 
target equilibrium states associated with the two environments can be introduced 
as L aj(t) ~ Lq{AJ{t)/J)- 1 /^ and L AT{r) ~ Lq(AT{t) /T^)- 1 ^ . The target equi- 
librium states are completely different at length scales much larger than the overlap 
lengths by the chaos effects. 5 ) During the continuous perturbations, the relaxational 
process which tends to grow the domains of W(t) and rejuvenation due to the chaos 
effects must be competing with each other. Naturally we expect that the system 
attains a certain stationary state characterized by a certain fixed size of the domains 
Lt(t*) with respect to the temporal working environment W(t). In other words 
the age of the system is pinned to a certain value r*. As proposed in Ref, 3 )' 8 ) such 
a length scale should play the role of a sort of renormalized overlap length which 
parameterize the heating/cooling rate effects in temperature-cycling experiments. 
Simplest possible scaling form for the stationary domain size is 

L t (t*) ~ L (AJ(t*)/J)- 1/<: L t (t*) ~ Lo(AT(T*)/T g )-W (1) 

It can be seen that by choosing slower perturbations p — > or vt — > 0, the pinned 
age increases r* — > oo and the stationary domain size increases Lt(t*) — > oo. 

It can be also seen that in the limit of slow perturbations p — > or v t — > 
actual change of the working environment becomes infinitesimal AJ(t*) — > or 
AT{t*) — ► 0. The later feature is quite unusual compared with analogues of more 
conventional cases such as systems under shear with shear rate 7. There one typically 
finds that relaxation times scales as r* oc j~ u with v < 1 so that actual deformation 
of the system 7 cx j 1 ~ u within the relaxation time r* can diverge as 7 — ► 0. 

During a continuous change of temperature the growth law Lx(t) would also 
change. However it is asymptotically irrelevant in the above scenario since Lt+at^) / L 
(T + AT/T) 1 ^ = O(l) . Thus for a given cooling/heating rate vt the relaxation 
time t* is defined at every local temperature T. The effective age r* of a system 
cooled down to a certain target temperature T m will be determined only by local 
cooling rate in the vicinity of T m and the cooling rates used at higher temperatures 
are irrelevant. This feature observed experimentally in spin-glasses 2 ) would hardly 
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be present in systems with activated dynamics but without the chaos effect. 
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Fig. 1. Spin autocorrelation function C(r + t w ,t w ) = (1/N) Yli=i (^i( T + £w)5i(tw)) under con- 
tinuous bond perturbation of strength p = 7.8125 x 10 -6 at T/J — 1.2. Here the averages 
over different realizations of thermal noises and sequences of bonds are denoted as (■ ■ ■ ) and 7-77 
respectively. The lines in the left figure are reference data of standard isothermal aging. 7 ' The 
system exhibits interrupted aging: waiting time t w dependence saturates. The tail parts of the 
autocorrelation functions are fitted to a simple exponential form j4exp(— t/r*(p)) by which the 
characteristic relaxation time r*(p) is extracted. In this example r* = 6013 (MCS) is obtained. 
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Fig. 2. The correlation Cbond(T*) between the set of bonds at two times separated by the re- 
laxation time r*(p) is shown in the left figure. It is defined as Cbond(i) = Jij(0)Jij(t)/J 2 = 
exp(— t/rbondip)) with the decorrelation time of the bonds given as Tbond(p) = — 1/ log(l — 2p) ~ 
1/p. As shown in the inset the relaxation time of the spins r*(p) grows with decreasing p but less 
rapidly than the decorrelation time of the bonds Tb on d(p)- In the right figure the scaling ansatz 
Eq. Qis tested. The open symbols are the reference data of Lt(t) (from Ref 7 '). The filled 
symbols are the overlap length a(AJ(r)/ J)" 1 ^ with AJ(t) = J^p e s(t) at the relaxation time 
t = T*{p) with p = 0.002 x 2~™ with n = 0, . . . , 7 (from left to right) obtained at T/J = 0.8, 1.2 
and 1.6. Here £ = 0.9 as obtained in Ref. 3 - 1 (see also 6 ') and chosen the prefactor a — 0.22. 



Numerical analysis- We performed MC simulations of the dynamics of the 
4-dimensional EA model (T g ~ 2.0J) under continuous bond changes starting from 
random initial configurations. In Fig. ^ we show the data of the spin autocorrelation 
function C(r + t w ,t w ). Apparently the system becomes stationary within a finite 
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time scale r* due to the continuous bond changes. The relaxation time r* extracted 
from the data are examined as shown in Fig. |2J It appears that the actual change 
of the bonds within the relaxation time r* of the bonds becomes negligible as the 
intensity of the driving p becomes smaller. As shown in the right figure of Fig. |2 
the relaxation time r* follows the expected scaling Eq. (^Q). 

Finally let us examine the fluctua- 
tion dissipation relation (FDR) as pro- 
posed by the dynamical mean field the- 
ory (DMFT). 9 ) As shown in Fig. El the 
fluctuation dissipation theorem is satis- 
fied at short time scales and violated at 
larger time scales in a non-trivial man- 
ner as predicted by the DMFT. Thus 
the stationary state can be considered 
as a glassy state. The result is very sim- 
ilar to that of the same model during 
isothermal aging 7 ) and driven dynamics 
induced by asymmetric coupling. 10 ) It 
must be noted however that the FDR 
systematically depends on p as those of 
isothermal aging systematically depends 
on t w so that the FDR in the asymptotic 
limit p — > may be different from the 

prediction of the DMFT. 7 ) , , , 

_ _ neat bath temperatures 1 . 

Discussion- To conclude we studied glassy stationary dynamics of spin-glasses 
under continuous temperature or bond changes. Let us consider a typical cooling rate 
vt = O(10~ 1 )(K/sec) used in " qunech" experiments, which is actually as slow as 
vt = 10~ 15 J/MCS since T g = O(10) K and the microscopic time scale is tq = 10~ 13 
(sec). The scaling ansatz Eq. (^Q) implies it yields r* of only order 1 — 10 (sec) at a 
target temperature T m = 0.7T g in the case of a canonical spin-glass AgMn sample, 3 ) 
which allows one to observe "clean" aging at T m . Indeed the effective age of the 
system after such slow cooling is found to be only O(10) (sec) experimentally. 3 ) 




Fig 



3. The ZFC susceptibility xzfc(t + 
t w ,tw) vs. the spin autocorrelation func- 
tion C(t + tw,tw) with tvt S> t* . In this 
example p — 2.5 x 10 -4 . Xzfc is measured 
as usual 7 ' with h/J — 0.1 but under con- 
tinuous bond changes. The lines represents 
FDT xzfc = (1-CQ/T associated with the 
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